Detection of valley polarization in graphene by a superconducting contact 
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Because the valleys in the band structure of graphene are related by time-reversal symmetry, 
electrons from one valley are reflected as holes from the other valley at the junction with a super- 
conductor. We show how this Andreev reflection can be used to detect the valley polarization of 
edge states produced by a magnetic field. In the absence of intervalley relaxation, the conductance 
Gns = (2e 2 /ft)(l — cosO) of the junction on the lowest quantum Hall plateau is entirely deter- 
mined by the angle between the valley isospins of the edge states approaching and leaving the 
superconductor. If the superconductor covers a single edge, = and no current can enter the 
superconductor. A measurement of Gns then determines the intervalley relaxation time. 

PACS numbers: 74.45.+C, 73.23.-b, 73.43.-f, 73.50. Jt 
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The quantized Hall conductance in graphene exhibits 
the half-integer quantization Gh = (n+ |0(4e 2 //i) char- 
acteristic of massless Dirac fermions 0, Q. The lowest 
plateau at 2e 2 /h extends to zero carrier density because 
there is no gap between conduction and valence bands, 
and it has only a twofold spin degeneracy because it lacks 
the valley degeneracy of the higher plateaus. The valley 
degeneracy of the lowest Landau level is removed at the 
edge of the carbon monolayer, where the current-carrying 
states at the Fermi level are located. Depending on the 
crystallographic orientation of the edge, the edge states 
may lie fully within a single valley, or they may be a lin- 
ear combination of states from both valleys [3, H| ■ The 
type of valley polarization remains hidden in the Hall 
conductance, which is insensitive to edge properties. 

Here we propose a method to detect the valley po- 
larization of quantum Hall edge states, using a super- 
conducting contact as a probe. In t he p ast, experimental 

HUSH and theoretical BESmillEl studies of the 
quantum Hall effect with superconducting contacts have 
been carried out in the context of semiconductor two- 
dimensional electron gases. The valley degree of freedom 
has not appeared in that context. In graphene, the ex- 
istence of two valleys related by time-reversal symmetry 
plays a key role in the process of Andreev reflection at the 
normal-superconducting (NS) interface A nonzero 

subgap current through the NS interface requires the con- 
version of an electron approaching in one valley into a 
hole leaving in the other valley. This is suppressed if the 
edge states at the Fermi level lie exclusively in a single 
valley, creating a sensitivity of the conductance of the NS 
interface to the valley polarization. 

Allowing for a general type of valley polarization, we 
calculate that the two-terminal conductance Gns (mea- 
sured between the superconductor and a normal-metal 
contact) is given by 



2e 2 

Gns = -t-(1 - cos 9), 
h 



(1) 



when the Hall conductance Gh = 2e 2 /h is on the lowest 
plateau [l5[ . Here cos 9 = ui ■ v-x is the cosine of the 
angle between the valley isospins V\ , of the states along 




FIG. 1: Three diagrams of a graphene sheet contacted by 
one normal-metal (N) and one superconducting (S) electrode. 
Edge states approaching and leaving the superconductor are 
indicated by arrows. The solid line represents an electron 
state (green: isospin v\\ blue: isospin 1/2), and the dashed 
line represents a hole state (red: isospin —v-i). 



the two graphene edges connected by the superconductor 
(see Fig. [TJ . If the superconductor covers a single edge 
(Fig.[TJi), then 9 = => Gns =0 — no current can enter 
into the superconductor without intervalley relaxation. If 
the superconductor connects different edges (Figs. Q]3,c) 
then Gns can vary from to 4e 2 //i — depending on the 
relative orientation of the valley isospins along the two 
edges. 

We start our analysis from the Dirac-Bogoliubov-De 
Gennes (DBdG) equation fil| 



H -n 
A* (i 



A 
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(2) 



with H the Dirac Hamiltonian, A the superconducting 
pair potential, and T the time reversal operator. The 
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excitation energy e is measured relative to the Fermi en- 
ergy /j,. Each of the four blocks in Eq. §2$ represents a 
4x4 matrix, acting on 2 sublattice and 2 valley degrees 
of freedom. The wave function ^ = (^ e ,^h) contains a 
pair of 4-dimensional vectors \& e and that represent, 
respectively, electron and hole excitations. 

The pair potential A is isotropic in both the sublattice 
and valley degrees of freedom. It is convenient to choose 
a "valley isotropic" basis such that the Hamiltonian H is 
isotropic in the valley degree of freedom 



H = 



(p + eA) ■ a 







(p + eA) ■ a 

vt <g) (p + eA) ■ a, 



(3) 



with v the Fermi velocity, p — (h/i)(d/dx,d/dy) the 
canonical momentum operator in the x-y plane of the 
graphene layer and A the vector potential corresponding 
to a perpendicular magnetic field B. The Pauli matrices 
<7j and Tj act on the sublattice and valley degree of free- 
dom, respectively (with <tq and tq representing the 2x2 
unit matrix). The time reversal operator in the valley 
isotropic basis reads 



T 





-ia,. 







C 



-(Ty ® <Jy)C, 



(4) 



with C the operator of complex conjugation. For later use 
we note that the particle current operator J — (J e , Jh) 
has electron and hole components 



J = v(t <g> <x, —t <g> a). 



(5) 



Substitution of Eqs. © and (g|) into Eq. © gives the 
DBdG equation in the valley isotropic form 



H+ - fi A 
A* fi-H- 



H± = vtq (g) (p ± eA) ■ a. 



(0) 
(7) 



We seek a solution in the normal region (where A = 0), 
at energies below the excitation gap A in the supercon- 
ductor. Electron and hole excitations cannot propagate 
into the superconductor at subgap energies, and the mag- 
netic field confines them in the normal region to within a 
magnetic length l m = yjTi/eB of the edge. We consider 
separately the edge states along the insulating edge of 
the graphene layer and along the interface with the su- 
perconductor. The edges are assumed to be smooth on 
the scale of l m , so that they may be treated locally as a 
straight line. We also assume that the superconducting 
coherence length £ = hv/Ao is small compared to l m , 
so that the effect of the magnetic field on the supercon- 
ductor can be neglected. 

The edge states at the insulating and superconducting 
boundaries are different because of the different bound- 
ary conditions. Using only the condition of particle cur- 
rent conservation, these have the general form [ItJ 



(8) 



with M. a unitary and Hermitian matrix that anticom- 
mutes with the particle current operator: 

M=M\ M 2 = l, M(n- J) + (n- J)M = 0. (9) 

The unit vector n lies in the x-y plane, perpendicular to 
the boundary and pointing outward. 

At the NS interface the matrix A4 is given by [l8[ 



M = 







M- 



NS 



M NS \ 
) 



M N s = t ® e 



i<p-\-ij3n-cr 



(10) 



with j3 = arccos(e/Ao) € (0, w) determined by the order 
parameter A = Aoe 1 ^ in the superconductor. 

The insulating (I) edge does not mix electrons and 
holes, so M. is block-diagonal with electron block M\ and 
hole block TM{F~ X . The boundary condition is deter- 
mined by confinement on the scale of the lattice constant 
i *C l m , so it should preserve time-reversal symmetry. 
This implies that Mi should commute with T. The most 
general matrix that also satisfies Eq. © is given by [l9[ 



M 



Mi 
Mi 



Mi 



(n±-tr), (11) 



parameterized by a pair of three-dimensional unit vectors 
v and n±. The vector n± should be orthogonal to n 
but v is not so constrained. Three common types of 
confinement are the zigzag edge, with v = ±z, n± = z; 
the armchair edge, with v ■ z = 0, n± ■ z = 0; and infinite 
mass confinement, with u = z, n± ■ z = 0. 

To determine the edge states we consider a local co- 
ordinate system such that the boundary is along the y- 
axis (so n = —x), and we choose a local gauge such 
that A = Bxy. The wave number q along the boundary 
is then a good quantum number. In order to simplify 
the notation we measure energies in units of hv/l m and 
lengths in units of l m . (Units will be reinstated in the fi- 
nal results.) Eigenstates of Eq. © that decay for x — > oo 
have the form 



9 ^=^[ciilZ-V' (12) 

*.({) = e-\? H(^)H ( , +s) V2-dO\ (13) 



-«(//- E)F (jU _ e)2/2 _i(0 



. (14) 



in the region x > (where A = 0). The function H a (x) 
is the Hermite function. The two-component spinors C e 
and Ch determine the valley isospin of the electron and 
hole components, respectively. 

The dispersion relation between energy e and momen- 
tum q follows by substitution of the state (|12[) into the 
boundary condition ([8]) . At the NS interface we take Eq. 
(fT0| for the boundary condition and obtain 

e[U +e (q)U- e (-q) + 1] 



U+s(q) - U-e(-q) 



Mi) 



H, 



(g) 

aH a 2 /2 -i{q) ' 



a 2 /2 



(15) 
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FIG. 2: Dispersion relation of edge states in graphene along 
the normal-superconducting interface, calculated from Eq. 
(|15[) for |e| <C Aq. The dotted lines are for /i = 0, the solid 
lines for p = 0Ahv/l m - 
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FIG. 4: Dispersion relation of states along the insulating edge, 
calculated from Eqs. (|16[1 and (|17[1 for /i = 0Ahv/l m and 
6 = 7r/2. The solid lines are the electron states (blue ej", red 
£e ), the dashed lines are the hole states (blue red e^). 
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FIG. 3: Cyclotron orbits of Andreev reflected electrons and 
holes. 



The solutions e n (q), numbered by a mode index n = 
0, ±1,±2, are plotted in Fig. O Notice that the 
dispersion relation has the inversion symmetry e(q) — 
—e(—q). Each mode has a twofold valley degeneracy, 
because the boundary condition (jlOjl is isotropic in the 
valley isospin v. The two degenerate eigenstates (labeled 
±) have Cf = c e \±u), = Ch\ ±f), with | ± v) eigen- 
states of V ■ T [21j . 

The expectation value v n — h~ 1 de n /dq of the veloc- 
ity along the boundary in the n-th mode is determined 
by the derivative of the dispersion relation. We see from 
Fig. [5] that the edge states all propagate in the same di- 
rection, dictated by the sign of B and (x. The velocity 
vanishes for \q\ — > oo, as the NS edge states evolve into 
the usual dispersionless Landau levels deep in the normal 
region. For q — > — oo the Landau levels contain electron 
excitations at energy e n = \/2(hv/l m ) sign — fi, 

while for q — > oo they contain hole excitations with 
£ n = \/2(hv/l m ) sign (n)yjn| + \i. For fi = the NS 
edge states have zero velocity at any q for |e| <C A . As 
illustrated in Fig.[3l the localization of the edge states as 
fi — > happens because for |e| > the electron and hole 
excitations move in opposite directions along the bound- 
ary, while for |e| < \fi\ they move in the same direction. 



Turning now to the insulating edge, we take the bound- 
ary condition (jTTJ) . For an edge along the y-axis we have 
n± = (OjSin^jCosf?). The valley degeneracy is broken 
in general, with different dispersion relations for the two 
eigenstates | ± v) of v ■ t. The dispersion relations for 
electrons and holes are related by e^(q) = —ef(—q). For 
sufficiently small /i there is one electron and one hole 
state at the Fermi level, of opposite isospins. (Note that 
electrons and holes from the same valley have opposite 
isospins.) We fix the sign of v such that | + u) is the 
electron eigenstate and | — v) the hole cigcnstatc. Wc 
find that £^{q) is determined by the equation 

f M+e {q) =tan(0/2), (16) 

while £~{q) is determined by 

/ M+e (g) = -cotan(e/2). (17) 

The dispersion relations plotted in Fig. 2] are for the case 
9 = tt/2 of an armchair edge. The case 9 = of a zigzag 
edge contains additional dispersionless states away from 
the Fermi level [H , but these play no role in the electrical 
conduction. 

To determine the conductance Gns we need to cal- 
culate the transmission matrix t of the edge states at 
the Fermi level. Edge states approach the superconduc- 
tor along the insulating edge L (with parameters i>i, 9\), 
then propagate along the NS interface, and finally return 
along the insulating edge I2 (with parameters 1/2, #2)- At 
sufficiently small /i each insulating edge I p supports only 
two propagating modes, one electron mode oc | + u p ) and 
one hole mode oc | — v p ). The NS interface also supports 
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two propagating modes at small p, of mixed electron- 
hole character and valley degenerate. The conductance 
is given by [13] 

2e 2 4e 2 

GnS = -T - (1 — T ee + The) = -r— Tft e , (18) 

with T ee = \t++\ 2 the probability that an electron in- 
cident along Ii returns along ±2 as an electron and 

The = \t h| 2 the probability that the electron returns 

as a hole. Because electrons and holes cannot enter into 
the superconductor, these two probabilities must add up 
to unity — hence the second equality in Eq. (| 18|) . (The 
factor of two accounts for the spin degeneracy.) 

Since the unidirectional motion of the edge states pre- 
vents reflections, the transmission matrix t from E to I2 
is the product of the transmission matrices t\ from Ii to 
NS and t% from NS to I2. Each of the matrices t p is a 
2x2 unitary matrix, diagonal in the basis | ± u p ): 

t p = + v p ){+v p \ + e<\ - v p )(-v p \. (19) 

The phase shifts (j) p , (j)' p need not be determined. Using 
I (I'll ± 1^2) 1 2 = |(1 ± Vi ■ i/ 2 ), we obtain from t — t 2 t\ the 
required transmission probabilities 

T he = l-T ee = ±(l-v 1 -v 2 ). (20) 

Substitution into Eq. (fl8|) gives our central result |T]). 

Referring to Fig.Q] we see that Gns = in the case (a) 
of a superconducting contact to a single edge (y\ = v 2 ) 
— regardless of whether the edge is zigzag or armchair. 
In the case (c) of a contact between a zigzag and an 
armchair edge we have V\ ■ v 2 = =>• Gns = 2e 2 //i. 



The case (b) of a contact between two opposite edges 
has i>»i = —v 2 => Gns = 4e 2 //i if both edges are zigzag; 
the same holds if both edges are armchair separated by 
a multiple of three hexagons (as in the figure); if the 
number of hexagons separating the two armchair edges is 
not a multiple of three, then 1/1-1/2 = 1/2 => Gns = e 2 /h. 

Intervalley relaxation at a rate T tends to equalize 
the populations of the two degenerate modes propa- 
gating along the NS interface. This becomes appre- 
ciable if TL/vq > 1, with L the length of the NS 
interface and vq = h~ 1 deo/dq ~ min(i>/2, \/2 pl m /h) 
the velocity along the interface. The density matrix 
p = /9 (1 - e~ TL / v °) + pie~ TL / v ° then contains a val- 
ley isotropic part p^ cx r with T ee = T e h = 1/2 and a 
nonequilibrium part po cx |i/i)(i/i| with T eei T e h given by 
Eq. ([20]) . The conductance then takes the form 

Gns = ^(l-e- TL 'v° cos 8). (21) 

A nonzero conductance when the supercurrent covers a 
single edge (9 = 0) is thus a direct measure of the inter- 
valley relaxation. 

In conclusion, we have shown that the valley structure 
of quantum Hall edge states in graphene, which remains 
hidden in the Hall conductance, can be extracted from 
the current that flows through a superconducting con- 
tact. Since such contacts have now been fabricated suc- 
ccsfully [12, [H[, we expect that this method to detect 
valley polarization can be tested in the near future. 
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